Draft version February 1, 2008 

Preprint typeset using LAT^X style emulatcapj v. 04/03/99 



DISCOVERY OF A MAGNETIC WHITE DWARF IN THE SYMBIOTIC BINARY Z 

ANDROMEDAE 

J. L. SOKOLOSKI AND LARS BlLDSTEN 
Department of Physics and Department of Astronomy, University of California, Berkeley, Berkeley, CA 
94720-3411; jeno@song.berkeley.edu, bildsten@fire.berkeley.edu 
Draft version February 1, 2008 

ABSTRACT 



OO 
ON 

On 
o 

Q 



We report the first result from our survey of rapid variability in symbiotic binaries: the discovery of a persistent 
oscillation at P — 1682.6 ±0.6 s in the optical emission from the prototype symbiotic, Z Andromedae. The oscillation was 
detected on all 8 occasions on which the source was observed over a timespan of nearly a year, making it the first such 
persistent periodic pulse found in a symbiotic binary. The amplitude was typically 2 — 5 mmag, and it was correlated 
with the optical brightness during a relatively small outburst of the system. The most natural explanation is that the 
oscillation arises from the rotation of an accreting magnetic (Bg J> 10 5 G) white dwarf. This discovery constrains the 
outburst mechanisms, since the oscillation emission region near the surface of the white dwarf was visible during the 
outburst. 

Subject headings: accretion — binaries: symbiotic — stars: magnetic fields — stars: oscillations — 
stars: rotation — white dwarfs 
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1. INTRODUCTION 

When the term "symbiotic star" was coined in the early 
1940s for the newly discovered pec uliar variable s tars with 
combination optical spectra (see Kcnyon 1986), Z An- 
dromedae was one of the prototypes. Today, it remains 



one of the mo3t frequently observed symbiotic systems 
(SS). The observations have rev ealed a complex eystom 



that is still not fully understood ( Mikolajewska fc Kcnyon 
1996). Most evidence indicates that the hot star in Z An- 
dromedae is a white dwarf (WD), and the work we present 
here supports that conclusion. This evidence includes ef- 
fective temperature estimates of the hot compone n t of ap 



proximately 10 5 K ( Fernandez-Castro et al. 1988) ; Mursct 



et al. 1991), an i nferred hot component radius of approx- 
imately .07.Bp) ( Fcrnandcz-Castro et al. 1988[ [Vlursct et 
al. 1991), an d a large radio nebula (Seaquist, Taylor, & 



Button 1984), which is not expected if mass transfer occurs 
via Roche lobe overflow onto a main s equence star. The 
binary has an orbital period of 759 days (Formiggini fc LcT| 
bowitz 1994| ; [Mikolajewska fc Kenyon 1996 ), and Schmid 
and Schild (1997) have used Raman line polarimetry to 
determine an orbital inclination of i 47 ± 12° and infer 
a mass for the hot component of 0.65 ± 0.28M Q (assuming 
a total system mass of between 1.3 and 2.3 M Q ) 



According to cur rent theories of binary evolution ( Yun- 
gelson et al. 1995), most white dwarfs found in symbi- 
otics should have evolved from stars with main sequence 
masses greater than about 1.5M . Highly magnetic WDs 
(Bs ;> 10 6 G, where B$ is the field at the stellar surface) 
appear to be preferentially for med by stars with main se- 
quen c e masses M ~ 2 — 4 M^ ( Angel, Borra, fc Landstrec^ 
1981; [gion et al. 1988 ), and so it is possible that the 



fraction of WDs that are magnetic is higher in symbiotics 
than in the field, where it is about 3 — 5% ( Chanmugam| 
1992). Given that there are at least 150 known SS, and 
that most of these contain white dwarfs, we expect some 
SS to contain white dwarfs that are magnetized at the 
level seen in DQ Hcrculis and AM Herculis cataclysmic 
variables (B s £ 10 5 G). Mikolajewski et al. (1990ab) and 



also Mikolajewski & Mikolajewska (1988) have invoked the 
presence of a magnetic white dwarf to explain the jets, 
flickering with possible QPO's, and large changes in the 
hot component luminosity in the symbiotic star CH Cygni, 
and this idea was later also adopted in the case of another 



Michalitsianos 



symbiotic, M WC 560 flTomov et al. 1992| 
et al. 1993| ). However, stable and repeatable oscillations 



like those detected in magnetic cataclysmic variables have 
until now not been seen in a symbiotic, and the prevalence 
of magnetic white dwarfs in SS is an important unknown. 

We are undertaking a long-term observational program 
to study the minute time scale photometric behavior of 
symbiotic binaries, expanding upon the work of Dobrzy- 
cka, Kenyon, & Milone (1996) and others. In §0, we 
present the first result from our survey, the discovery of a 
28-minute oscillation in the B band emission from Z An- 
dromedae. This was the only strong oscillation found in 
the preliminary analysis of 20 objects. Results from the 
complete survey will be presented in a future paper. In 
we interpret this oscillation as due to accretion onto a 
magnetic, rotating white dwarf. The fact that the oscilla- 
tion was detected during a recent outburst, as well as once 
the source had returned to quiescence, has implications for 
outburst models, as we outline in our conclusions (§^). 

2. OBSERVATIONS AND RESULTS 

We observed Z Andromedae on seven occasions sepa- 
rated by two to four weeks each from 1997 July to 1997 
December, and then once again in 1998 June, with the 
1-meter Nickel telescope at UCO/Lick Observatory. The 
observations ranged in length from approximately 4 hours 
on a single night in 1997 July, to approximately 7 hours 
per night three nights in a row in 1997 November, for a 
total of 76 hours of observing on 13 nights spanning 1 year 
(see Table 1). The 2048 x 2048 pixel, 6.3 x 6.3 arcmin, un- 
thinned LORAL CCD currently in Lick's dewar #2, and 
a Johnson B filter were always used. The first observation 
fortuitously occurred about 1 month after the peak of an 
optical outburst (AV w 1), and the subsequent observa- 
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tions in 1997 took place as the optical flux declined. At 
the time of the 1998 June observation, the optical flux had 
returned to its pre-outburst quiescent value. A 2.2 year V 
band light curve is shown in Figure 1, with the times of 
our observations marked. 
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Fig. 1. — Long term V band light curve of Z *nfli-p' r '° rla °, fr "p 
the American Association of Variable Star Observers (Mattci 1998). 
The times of our 8 observations are marked with arrows; sec Tabic 
1 for the corresponding dates. 

At the time of the first observation, on 1997 July 8, the 
binary was oriented with the WD in front of the red giant, 
from the observer's perspective (i.e. the orbital phase of 
the binary was 0.5, where phase 0.0 corresponds to photo- 
metric minimum in quiescence and spectroscopic conjunc- 
tion). By 1997 December 2, the binary had moved through 
almost one quarter of its orbit to phase 0.7, where both 
stars are roughly equidistant from the observer, and in 
1998 June, the WD was behind the red giant. 

Data reduction was performed using IDL software based 
on standard IRAF routines. Source counts for each image 
were extracted from a circular aperture with a radius of 
8 to 14 arcsec, and the background was estimated from 
a surrounding annulus. For each light curve, the extrac- 
tion region was chosen to be much larger than the seeing 
so that any variability due to source counts falling out- 
side the extraction region (as the seeing or guiding qual- 
ity changed) was small compared to systematic errors. Z 
Andromedae is bright enough that even with large extrac- 
tion regions, the Poisson errors from sky background are 
usually not significant. Several representative light curves 
from our observations are shown, in chronological order, 
in Figure 2. There is one other bright star in the field of Z 
Andromedae (at J2000 coordinates 23 33 24, +48 45 38), 
but it is also variable, so it was not used as a compari- 
son star for differential photometry (except for one night 
in October when its amplitude of variability was low, and 
one night in November when thick clouds were present). 
Therefore, although every attempt was made to perform 
observations on clear nights, some observations were af- 
fected by high clouds. Data points affected by radiation 
events ( "cosmic rays" ) were removed when they could be 
identified, and the light curves were corrected for atmo- 



spheric extinction. In addition, we divided most of the 
light curves by a 3rd order polynomial in order to remove 
residual atmospheric effects (<j 1% effect). Note that this 
may have also removed any intrinsic variability on a time 
scale comparable to the length of the observation. This 
polynomial fitting was not performed for the Aug. 2 and 
Aug. 30 data because of the presence of flare-like variabil- 
ity in the light curves (see Figure 2). 
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Fig. 2. — Four B Band light curves of Z Andromedae. The light 
curves have been corrected for atmospheric extinction, and normal- 
ized to a typical value of the count rate. The 28-minute oscillation 
is clearly visible in at least 3 of the 4 light curves. The additional 
variability on August 2 and August 30 could be intrinsic to Z An- 
dromedae, but the lack of a constant comparison star prohibits us 
from ruling out atmospheric origins. 

Power spectra corresponding to the light curves in Fig- 
ure 2 are shown in Figure 3. The most striking, persistent 
feature is the peak at 0.6 mHz, corresponding to a period 
of 28 minutes. A smaller, but still significant peak is also 
present at twice this frequency in the power spectra of the 
July 8 (not shown), August 2, August 30, and September 
13 light curves. The 28-minute oscillation was significantly 
detected in all 8 observations. The power spectrum of the 
other bright star in the field does not show the feature 
at 0.6 mHz, confirming that the oscillation detected in Z 
Andromedae was neither instrumental nor atmospheric in 
origin. 
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Fig. 3. — Power spectra for the light curves in Figure 2. The 
power is plotted in units of mean high frequency power (which is 
"white"), and the feature at 0.6 mHz corresponds to an oscillation 
period of 28 minutes. A smaller feature at 1.2 mHz is also seen 
(and significantly detected) in the top three panels. The features at 
frequencies lower than 0.6 mHz are not repeated from one observa- 
tion to another, and could be due to atmospheric changes or source 
variability. 

2.1. Timing Analysis 

Although the feature at 0.6 mHz is detected in the power 
spectra of Z Andromedae, the precise value of the oscilla- 
tion period and its uncertainty is best determined in the 
time domain. Visual inspection of the light curves sug- 
gests that the signal is not a simple sinusoid, and the de- 
tection of a harmonic in several of the power spectra con- 
firms this impression. By using time domain epoch folding 
techniques, we need not make any assumptions about the 
shape of the pulse profile, and all the signal "power" will 
be located at a single period. 

To perform this analysis, we used the phase disper- 
sion minimization (PDM) technique originally described 
by Stellingwerf (1978). The PDM technique consists of 
folding a light curve at a range of periods and computing 
the mean pulse profile, and the scatter of the data points 
about this profile, for each period. Each data point is as- 
signed a phase <fi = t mod P, where t is the time of the 
measurement from some initial time and P is the fold pe- 
riod, and binned accordingly. In our case the number of 
phase bins ranged from 10 to 20, and was chosen to be as 
large as possible while still ensuring that each phase bin 
contained enough data for our statistics to be valid (at 
least 10 data points). With a large number of points in 
each bin, the standard PDM statistic is just x 2 , 
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where n& is the number of phase bins, n, is the number of 
points in bin i, Xij is the jth point in the zth phase bin, 
TOi = n^ 1 ^2 is the mean of all points in bin i, and cr^ is 
the uncertainty on . In our treatment the <Ty were usu- 
ally dominated by the Poisson errors on the source counts. 
The mean pulse profile will be rather flat for fold periods 
far from the true period, and the points in each phase bin 
will have a large variance, causing \ 2 to be large. For a 
fold period close to the true period, the mean profile will 
approach the true pulse profile, the variance of the points 
in each bin will be small, and \ 2 will decrease. The best 
estimate of the true period is found by minimizing % 2 . 

For a data stream with Gaussian noise and a superim- 
posed oscillation, Xmim the minimum value of x 2 , should 
be approximately equal to the number of degrees of free- 
dom (in this case N — rib, where N is the total number of 
points). In other words, the reduced \ 2 is approximately 
equal to 1, indicating that the fit of the data to the mean 
profile is good. If the errors are normally distributed, there 
is also a simple relationship between A% 2 above the mini- 
mum and the level of confidence that the true period lies 
within the range that produce x 2 < (Xmin + ^X 2 )- For 
models with only one free parameter, like ours, A\ 2 = 1 
corresponds to a 68.3% confidence level, and A\ 2 = 2.71 
corresponds to a 90% confidence level. If underlying red 
noise is present at the period of interest, the task of identi- 
fying an oscillation and measuring its parameters is much 
more difficult. Very long baselines or repeated observa- 
tions are then needed to characterize the underlying vari- 
ability. For Z Andromedae, any red noise is at longer 
periods than the detected oscillation, so we did not have 
this added complication. 

For data with noise properties that are not precisely 
Gaussian, the PDM technique can still be used, with the 
help of simulations to determine the correct relationship 
between A\ 2 above the minimum and confidence levels 



(van der Klis, private communication; Press et al. 1992). 
We performed such simulations for each light curve by re- 
peatedly injecting a fake signal with a known period into 
the data, and then examining the distribution of periods 
resulting from the PDM method. 

The period of the oscillation as determined from each 
observation individually is shown in the last column of 
Table 1, where the quoted errors are roughly 68% confi- 
dence limits. The measurements from observations with 
more than one night of observing are more precise than 
single night observations because of the longer baseline. 
The period measurements from all observations are con- 
sistent, indicating that the period was stable to within 
less than 15 seconds, or 1% of the period, for 1 year and 
during an outburst of the system. More accurate deter- 
mination of the oscillation period in Z Andromedae by 
connecting the data from adjacent observations will allow 
for important orbital time delay measurements. Given the 
system inclination of i = 47° and taking the total system 
mass to be M to t — 2Afp, w hich is typical for a symbi- 



otic (Schmid & Schild 1997, and references therein), the 
light travel time across the WD orbit is iflwD sin i)/c » 
12.2min (sini/0.73) (Af tot /2M Q ) 1 / 3 (1 + M WD /M RG )-\ 
where awD is the distance from the WD to the center 
of mass, Mwd is the mass of the WD, and Meg is the 
mass of the red giant. 
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Table 1 
Observation Log and Results 



Observation 


Date, U.T. 


MJD 


Orbital 


Obs. Length 




At 


Count Rate 


Oscillation 


Number 




-2450000 


Phasc a 


(hrs) 


(sec) 


(sec) 


(xl0 4 c/s) 


Period (sec) 


1 


1997 Jul 8 


637 


0.50 


3.7 


18 


40 


2.9 


1676 ± 19 


2 


1997 Aug 2 


663 


0.54 


0.0 


30 


58 


2.3 


1686 ± 11 


3 


1997 Aug 30 


691 


0.58 


6.5 


30 


57 


1.9 


1695 ± 15 


4 


1997 Sep 13 


705 


0.59 


8.0 


38 


63 


1.9 


1682.2 ± 0.7 




1997 Sep 14 


706 


0.60 


5.1 


65 


90 






5 


1997 Oct 5 


727 


0.62 


6.5 


50 


78 


1.6 


1684.2 ± 1.5 




1997 Oct 6 


728 


0.62 


6.2 


52 


79 






6 


1997 Nov 1 


754 


0.66 


6.2 


40 


68 


1.4 


1682.7 ± 1.0 




1997 Nov 2 


755 


0.66 


7.7 


65 


93 








1997 Nov 3 


756 


0.66 


7.1 


40 


68 






7 


1997 Dec 2 


785 


0.70 


5.8 


110 


138 


1.2 


1679 ± 16 


8 


1998 Jun 28 


993 


0.97 


3.4 


200 


228 


0.5 


1682.0 ± 3 




1998 Jun 29 


994 


0.98 


3.7 


200 


228 







a From the Formiggini & Leibowitz (1994) ephemeris: Min(vis) = JD2442666(±10) + 758.8(±2)S, where E is the number 
of orbital cycles. 



Note. — t exp is the integration time, and Af is the time between integration starts, which is equal to the integration time 
plus the CCD readout and processing time. 



The peaks at twice the fundamental frequency in the 
early data indicate that the pulse profile deviates from a 
sinusoid. Pulse profiles created by folding the light curves 
from observations #3, #4, and #7 at 1683 seconds are 
shown are Figure 4. The pulse fraction decreased mono- 
tonically as the outburst decayed, from ps 5 mmag peak- 
to-peak in 1997 July and August (observations #1 - #3) 
to k 2 mmag peak-to-peak in 1997 December (observation 
#7). In 1998 June the oscillation was detected at the 2 
mmag level. 
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Fig. 4. — Z Andromedae pulse profiles at P = 1683 s from 1997 
July 7, 1997 August 30, and 1998 June 28-29. Error bars are left off 
for clarity of viewing, but 1 sigma is about 0.0003 in the normalized 
units plotted, as shown by the mark in the center of the plot. The 
data are repeated for clarity. The decrease in pulse amplitude for 
later observations is significant. 

3. THE CASE FOR MAGNETIC ACCRETION 

We interpret the 28-minute oscillation as the result of 
rotation of a white dwarf which has a strong enough mag- 
netic field to channel the accretion flow onto i ts mag- 



netic pola r caps, as in the DQ Herculis systems (Patter- 
son 1994). Non-radial g-mode pulsations of a hot WD 
that is similar to a planetary nebula nucleus (PNN) is an- 
other plausible explanation of the oscillation, especially 
since g-mode pulsations with periods close to 28 minutes 
have been observed in several PNN (Ciardullo & Bond 
1996). However, these systems are multiperiodic and have 
frequencies that change on month long time scales. The 
Z Andromedae emission oscillated at only a single, con- 
stant frequency for an entire year, as well as throughout 
an outburst during which conditions in the WD envelope 
presumably changed significantly. Therefore, we conclude 
that WD g-mode pulsations are unlikely to be the cause 
of the oscillation. The period of the oscillation is too long 
to be due to an acoustic (p-mode) pulsation in a WD, and 
too short to be due to a g-mode pulsation in a main se- 
quence star with M w 0.65M Q . A p-mode pulsation in 
a main sequence star is not formally ruled out, but again 
one would expect more than a single mode to be present. 
Therefore, the period, its stability and coherence, and the 
fact that only one period is detected all support the WD 
magnetic dipole rotator model. 

The minimum magnetic field strength at the dipolar cap, 
Bs , that is needed to funnel the accretion onto the star can 
be roughly estimated by requiring that the magnetospheric 
radius, r mag , be larger than the white dwarf radius, R. At 
the magnetospheric radius, the magnetic field pressure is 
comparable to to the ram pressure of the in-falling ma- 
terial, giving the standard r mag (fi 4 /2GMwdM 2 ) 1 ^ 7 , 
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wher e p. = Bc,R 3 /2 is the magnetic dipolc. This then outburst (|Fcriiaiidcz-Castro et al. 1995 , Mikolajewska & 



leads ]Q a minimum field 



Kenyon 1996 J! Mikolajewska & Kenyon (1996) deduced 



B s ^3x 10 4 G 



10 cm 
R 



M 



10- 8 M© yr- 1 



M 



WD 



1/2 



1/4 



(2) 



0.65M, 



If the WD has been spun up so that it is in a rotational 
equilibrium, with its spin period, P s , equal to the Kepler 
period at r mag , then the magnctosphcric radius would be 



( GMwd P* 

r ' nag ~ ^ 4tT 2 



2\ V3 



(3) 



that the radius of the hot component increased by a fac- 
tor of ~ 100 during previous outbursts of Z Andromedae. 
However, they also noted a few problems with the shell 
flash/photospheric expansion model. The Hell emission 
lines evolve in a different manner than other emission lines 
during outbursts. Therefore, the outburst spectra are in- 
consistent with an evolving single temperature model for 
the WD. Another problem for thermonuclear runaway and 
steady burning shell expansion models is the time scales. 
It is difficult to reconcile theoretical photospheric expan- 
sion time scales a nd shell flash recurrence time sc ales with 
the observations ( Mikolajewska fc Kenyon 1992| , and ref- 



0.26.R 







M WD \ 1/3 ( I'.. 



0.65M 



28 min 



2/3 



and the field strength needed to have the magnetosphere 
at this radius is roughly 



B s « 6 x 10 6 G 



M 



1/2 



10- 8 M Q y^ 1 



10 cm 



R 



(4) 



M 



WD 



5/6 



v 0.65M Qy 

The time it takes for the white dwarf to reach this 
rotational equilibrium is t S pin-up = 2irI/NP s w 5 x 
10 4 yr (28 min/P s )(10- 8 M s yr^/M), where I « 
Mwp R 2 /5 is the WD moment of inertia and N 



M(G kwm 



a/2 



erences therein), especially for a low-mass WD (although 
see Sion & Ready [1992]). Our detection of an oscillation 
from a region that would be hidden by an expanded photo- 
sphere is another phenomenon that is difficult to reconcile 
with models involving photospheric expansion. 

Our observations do not provide information about the 
temperature evolution of the hot component, so it is pos- 
sible that the 1997 outburst was significantly different 
from previous outbursts. The 1997 outburst was smaller 
and more asymmetric than either of two well-studied out- 
bursts in 1984 and 1986, which rose to V « 9.6 and 
V s» 9.1 respectively, compared to V « 9.7 for the 
1997 outburst. Based upon multiwavelength observations, 
Fernandez-Castro et al. (1995) suggested that the 1984 
and 1986 events were similar, but that a less massive shell 
was ejected during the 1984 outburst. The 1984 event pro- 
duced a smaller increase in opacity (Fernandez-Castro et 



is the accretion torque. This spin- up al. 199J ), so a correlation between V at the outburst peak 



time is shorter than the lifetime of the red giant, so it is 
likely that the system has reached this form of equilibrium. 

4. CONCLUSIONS AND IMPLICATIONS FOR THE 
OUTBURST MECHANISM 



and the nature of the outburst might exist 

Another outburst mechanism that has been discussed 
for SS, although usually not for systems th at contain WD s, 
is thermal accretion-disk instabilities (DI; Duschl 1986a b) 



Ou: survey has so far yielded one persistent periodic os- 



likc t hose that lead to dwarf novae eruptions in CVs ( Dsafc 



cillator. The oscillation was detected on all 8 occasions 
when the source was observed over the course of one year, 
and the period, P = 1682.6 ± 0.6 s, was stable to within 
our measurement errors. We interpret this oscillation in 
terms of magnetic accretion onto a rotating WD. This de- 



1996). DI models have not been considered prime candi- 



tectio i is the first of its kind for a symbiotic, and it comes 



dates for explaining the outbursts in WD SS for several 
reasons. First of all, there is little direct evidence for disks 
around WDs in SS. Disks are not needed in spectral fits 
(Murset et al. 1991) and double peaked line profiles can- 
not be d efinitively linked to disk emission (Robinson et 
al. 1994). Furthermore, disk instabilities alone may not 



from an object, Z Andromedae, in which no other phenom- 
ena thought to be associated with magnetism have been 
observed. Outburst mechanisms need to be reconsidered 
in light of this discovery, and as we now elaborate, accre- 
tion disk instabilities look to be a promising source for the 
outbursts. 

The detection of an oscillation that originates at the 
WD surface during an outburst has serious consequences 
for models of the outburst mechanism in Z Andromedae. 



Most models ( Mikolajewska fc Kenyon 1992| , and refer- 
ences therein) invoke dramatic expansion of the WD pho- 
tosphere, for example as the result of a thermonuclear 
shell flash or a change in M onto a nearly stably burn- 
ing hydrogen layer. Evidence for such expansion and the 
subsequent decrease in the effective temperature of the 
WD includes a decrease in the strength of high ioniza- 
tion state emission lines, line broadening, increased opac- 
ity as measured by line ratios, the appearance of an A-F- 
type spectrum, and direct luminosity estimates, all during 



provide sufficient ener gy to explain the observed flux in- 
creases (Kenyon 1986). Disk instability models can, how- 
ever, produce time scales that are more in accordance with 
the durations and recurrence times seen in SS outbursts 
than thermonuclear runaway models. We will explore the 
possibility that disk instabilities play an important role in 
SS outbursts more fully in a separate paper. 

There are several points that are important to note 
here, however. Most importantly, if a large disk does ex- 
ist around the WD in Z Andromedae, it would be ther- 



mally unstable (Duschl 1986b; Meyer-Hofmeister 1992). 
Secondly, during a Dl-induced outburst, the emission re- 
gion close to the WD could remain exposed, as appears 
to be the case during the most recent outburst of Z An- 
dromedae. Finally, the presence of a quasi-steady burning 
layer on the WD may affect the energetics of the outburst 
resulting from a disk instability. 
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